Long range Casimir force induced by transverse electromagnetic modes 
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We consider the interaction of two perfectly conducting plates of arbitrary shape that are inside 
a non-simply connected cylinder with transverse section of the same shape. We show that the 
existence of transverse electromagnetic (TEM) modes produces a Casimir force that decays only as 
1/a 2 , where a is the distance between plates. The TEM force does not depend on the area of the 
plates and dominates at large distances over the force produced by the transverse electric (TE) and 
transverse magnetic (TM) modes, providing in this way a physical realization of the 1 + 1 dimensional 
Casimir effect. For the particular case of a coaxial circular cylindrical cavity, we compute the TE, 
TM and TEM contributions to the force, and find the critical distance for which the TEM modes 
dominate. 



I. INTRODUCTION 



In the last years, there has been an increasing interest in the Casimir effect The new generation of experiments 
allowed a precise determination of the Casimir force, and stimulated theoretical calculations of the forces for 
'■ different geometries, including finite temperature and conductivity corrections. 

The analysis of the dependence of the Casimir force with the geometry is therefore both of theoretical and experi- 
mental relevance. In this paper, we will point out that, in non-simply connected electromagnetic cavities, the presence 
of TEM modes produces an additional contribution to the Casimir force, that it is independent of the section of the 
cavity and decays slower than the contributions of TE and TM modes. As far as we know, this is the first example 
' in the literature that illustrates the relevance of the TEM modes in the static Casimir effect (for its relevance in the 
, dynamical Casimir effect see )■ 

• Concretely, we will consider the interaction between two identical perfectly conducting plates that are inside a 
very long cylinder of the same section, that is also perfectly conducting. Thesegeometries are usually referred to as 
7—1 ■ "Casimir pistons" Q, and have received considerable attention recently @, 0, [J One of the reasons that triggered 
these investigations was the reconsideration of the Casimir energy for rectangular boxes, since repulsive forces have 
O^N been predicted when considering only the zero-point energy of the internal modes of the rectangular cavities @ . The 
validity of these results has been disputed [j| for at least two related reasons: the omission of the contribution of the 
^ , external modes, and the ambiguity in the renormalization of the divergent quantities (however, there is no consensus 
in the literature on these issues, see for instance 10]). In any case, the advantage of the pistons is that, as long as 
the surfaces are perfectly conducting, one can compute the Casimir energy and forces unambiguously and without 
considering the external modes to the cavity. The new aspect that we will introduce in this paper is the consideration 
of non-simply connected cavities (see FigflJ, allowing for the existence of TEM modes. 

The paper is organized as follows. In Section II we will describe the different contributions to the zero-point energy 
in a non-simply connected cavity . We will see that the Casimir energy for TEM modes is equivalent to that of a 
massless scalar field living in 1 + 1 dimensions and satisfying Dirichlet boundary conditions on the plates. Moreover, 
the Casimir energy for TE and TM modes is equivalent to that of a set of massive, 1 -I- 1-dimensional scalar fields, 
with the masses determined by the eigenfrequencies associated to the transverse section of the cavity. We will discuss 
the behaviour of the force when the distance between the plates is much larger or much smaller than the transverse 
dimensions of the cavity, and conclude that the TEM force dominates above a critical distance. In Section III we 
will present a detailed analysis of the particular case of a coaxial cavity of circular section. We will evaluate the 
contributions of TE and TM modes to the Casimir energy for this geometry using a combination of the analytical 
result for a 1 + 1 massive field and Cauchy's theorem to perform the summation over the effective masses. We will 
also compare these contributions with that coming from TEM modes, and find the critical distance as a function of 
the radii of the inner and outer cylindrical shells. Section IV contains our final remarks. 



II. CASIMIR ENERGY IN NON-SIMPLY CONNECTED CAVITIES 



Let us consider a very long electromagnetic cylindrical cavity, with an arbitrary section. We will assume that the 
cavity is non-simply connected, i.e. that there is a second cylinder, also of arbitrary section, inside the larger one (see 
FigfT]). The cavity is the annular region between the two cylinders and contains two plates (pistons) separated by a 
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distance a (the pistons cover only the annular region between the cylinders). All surfaces are perfectly conducting. 
The z-direction is the axis of the cavity, and we will denote by xj_ the coordinates in the transverse sections. 




FIG. 1: Two pistons separated by a distance a inside a non-simply connected cavity (the annular region between two cylinders 
of arbitrary section). The system is enclosed between another two plates separated by a distance L 3> a. All surfaces are 
perfect conductors. 

At the classical level, the electromagnetic field admits a description in terms of independent TE, TM and TEM 
modes, which are defined with respect to z-dircction. This is possible due to the particular geometries we are 
considering, that have an invariant section along the z axis. The TE and TM electromagnetic degrees of freedom can 
be written in terms of two different vector potentials Ate and Atm with null divergence and z-component ■ The 
TE electric and magnetic fields are given by 

Ete = —Ate ; Bte = V x Ate, (1) 
while the TM fields are given by the dual relations 

Btm = Atm ; Etm = V x Atm- (2) 
The vector potentials can be written in terms of the so called (scalar) Hertz potentials [H, [l2j as 

Ate = z x V^te, 

A T m = z x V0 T m, (3) 

The Hertz potentials </>te and 4>tm satisfy the wave equation with Neumann and Dirichlet boundary conditions on 
the lateral surfaces, respectively, and the opposite boundary conditions on the pistons. The eigenfunctions can be 
chosen of the form 



and 



<hE(t,z,x±) = e-™'sin( Vte(x ± ) (4) 

a 

TUT Z 

<pTu(t, z, xj_) = e~ lwt cos( )<Ptm(x_l) , (5) 

a 

where n is a non-negative integer and </?te,tm are eigenfunctions of the transverse Laplacian 

V5_(/?TE,TM = — A 2 (y9TE,TM ■ (6) 

Therefore, the eigenfrequencies associated to the TE and TM modes are 



3 



where Xf, N and xi D are the eigenvalues of Eq[S] when the eigenfunctions satisfy Neumann and Dirichlet boundary 
conditions, respectively. 

When the cylindrical cavity is non-simply connected, in addition to the TE and TM modes one should also consider 
the TEM modes, for which both the electric and magnetic fields have vanishing z component. Working with the usual 
vector potential A, the TEM solutions are of the form 

A(x ± ,z,t) = A ± (ic ± )<p TEM (z,t), 
E = -(c\0tem) A ± , 

B = (<9 2 </>tem) z x Aj_. (8) 

where </>tem(z, t) is an additional scalar field. The transverse vector potential has vanishing rotor and divergence, 
and zero tangential component on the transverse surfaces. Therefore, A_l is a solution of an electrostatic problem in 
the two transverse dimensions (in hollow cylindrical cavities the transverse potential vanishes and TEM modes do not 
exist). The scalar field </>tem satisfies Dirichlet boundary conditions on the longitudinal boundaries z = and z — a, 
and the longitudinal wave equation (df — <9^)0tem = 0- Thus, the eigenfrequencies of the TEM modes are 

w™ = nw/a . (9) 

In order to obtain the Casimir energy we introduce the regularized quantities 

E Icg (a) --^^ e " + - 2_^\ W k,n e k -" + W k,n e k - n ) = ^rcg W + ^rcg W + ^reg W ' i 1U ) 
n n,k 

and two additional pistons separated at a very large distance L, enclosing the system. The physical Casimir energy 
is defined as the difference 

E(a) = E Teg (a) + 2E Ieg (^^) - 3£ reg (|) (11) 

in the limit when the cutoff a tends to zero. Note that one can compute independently the TE, TM and TEM 
contributions to the energy. Note also that, as the pistons only cover the annular region between the cylinders, the 
internal modes of the smaller cylinder will be irrelevant for the interaction between plates. 

To proceed, we note that the Casimir energy for this geometry is formally equivalent to that of a set of scalar fields 
living in 1 + 1 dimensions and satisfying Dirichlet boundary conditions at z = and z — a. Indeed, Eq[3] implies that 
the TEM Casimir energy is equivalent to that of a massless scalar field, and the result is very well known [l3j 

E TEM {a) = _JL_ (12) 

12a 

From Eql3 we see that the TE and TM Casimir energies correspond to that of a set of massive scalar fields, with 
masses given by XkN & n d Afc n. T he Casimir energy E m for a field of mass m in 1 + 1 dimensions has been computed 
previously by many authors [lij usin g di fferent methods of regularization (see in particular Ref. Q for a calculation 



with an exponential cutoff). It reads [1 



1 ^ mK^lma) 

E m {a) = -—^ , (13) 

i=i 



where K\ is the modified Bessel function of the second kind. Using this result and the analogy between the TE and 
TM eigenfrequencies (Eq.[7]) with the eigenfrequencies of massive scalar fields in 1 + 1 dimensions, we can easily obtain 
the TE and TM contributions to the Casimir energy in the cylindrical cavity 



+ oo 



E TE 



(a) + E™(a) = — ^ ( X! ^ kNKl ^^ kN °^ + ^2 ^kpKi(2l\ kD a) \ 



2tt ^ \ ^ I I 

i=i 



This equation has been previously obtained in Ref. Q using a different method. We stress that the formula is valid 
for a cavity of arbitrary section. The energy can in principle be computed through a numerical evaluation of the 
eigenvalues of the transverse Laplacian. Alternatively, as we will describe in the next section, the summation over 
the eigenvalues can be performed using Cauchy's theorem. The force between pistons is easily obtained taking the 
derivative of the energy with respect to a. 
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Let us now discuss some generic properties of the different contributions to the Casimir energy. At small distances, 
when the separation between pistons is much smaller that the transverse dimensions of the cavity, one expects the 
proximity force approximation (PFA) to describe accurately the contributions of TE and TM modes. Indeed, using 
heat kernel techniques it can be shown 0] that, in this limit 



E TE (a)+E™(a)*-^, (15) 



A 

720 a? 

where A is the area of the transverse sections. This is of course the well known result for parallel plates. It is worth 
to stress that the geometric properties of the transverse section, as the area, are contained in the eigenvalues XkD and 
AkN, which play the role of the masses of the fields in the 1 + 1 dimensional analogy. 

On the other hand, in the opposite limit we have Xa ^> 1. The TE and TM contributions to the Casimir energy are 
dominated by the lowest eigenvalue Xmin, and have the typical exponential suppression associated to massive fields, 
i.e. 



E TE (a) + E™ (a) a -J b^L e -^ IN a (16) 

where n is the multiplicity of the eigenvalue. 

From these results, we conclude that, at small distances, the total Casimir energy is dominated by the TE and TM 
contributions: it behaves as 1/a 3 as for parallel plates, and it is proportional to the area of the pistons. For distances 
larger than a critical value, a > a c , the TEM is the leading contribution, and gives a long range Casimir energy that 
decays only as 1/a. This contribution, typical of a massless scalar field in 1 + 1 dimensions, is non-extensive, i.e. does 
not depend on the area of the pistons. The value of a c depends of course of the particular form of the transverse 
section. 

There are some additional properties which can be obtained using dimensional analysis. Let us denote by l\ and 
I2 the typical lengths associated to the sections of the internal and external cylinders of the cavity, respectively. On 
dimensional grounds we expect 

E TE +E TM = l f (ah\ _ 



If the critical distance is defined by 



then wc have that 



a c l\\ 7r 
h'h) = ~12 



(18) 



a c = l 2 g(h/l2). (19) 

In the particular case 1% <C a c , I2, on physical grounds we expect the functions / and g defined above to have well 
defined limits: 

£T£ + E ™ K M^°) ' ac = ^(0). (20) 



Indeed, this limit can be achieved by inserting a thin wire inside a hollow cavity, so that the TE and TM modes of the 
hollow cavity are not disturbed, and so the TE and TM contributions to the zero-point energy are almost independent 
of the presence of the wire. In this situation the critical distance becomes a linear function of l 2 . We will confirm this 
property in the particular example described in the next section. 

The existence of the long range Casimir TEM force is of conceptual interest. One can wonder whether is it also 
relevant from an experimental perspective, i.e. if there is a chance of measuring this force in future experiments. 
There are two major limitations: on the one hand, being non-extensive in the area of the plates, the absolute value of 
the force is very small, and therefore it could only be measured at extremely short distances. Moreover, in this regime, 
the Casimir force would be dominated by the TE and TM contributions, unless the area of the plates is also small. 
We therefore address the question of whether the TEM contribution to the force can be a significant fraction (a) of 
the total Casimir force in the PFA regime, i.e. we are interested in a configuration in which F TEM > a(F TE + F™). 
For the sake of concreteness we study the case where the geometry of the plates is such that its area may be written 
as 6(^2 — '1), where b is some coefficient (this occurs, for instance, with a circle and any regular polygon). Within 
these assumptions, the conditions imposed to the forces become the following conditions in the parameters, 

A 20 

&« — <—, 21 

or ira 



■5 



where the first inequality follows from the validity of the PFA. This shows that, in principle, there is a region in 
the parameter space where the TEM force could be a significant fraction of the force, even in the proximity limit. 
However, we stress again that the smallness of the TEM force implies that its measurement is presently extremely 
difficult. 



III. COAXIAL CYLINDRICAL CAVITY OF CIRCULAR SECTION 



In this Section we present and solve an explicit example of a non-simply connected cavity. We compute the Casimir 
energy for a cavity formed by two concentric perfectly conducting circular cylinders closed by two (also non-simply 
connected) plates in its extremes. The configuration is similar to that in Fig[T]but with both inner and outer cylinders 
having circular sections of radii r\ and r 2 , respectively. We also study the Casimir force between the plates. 

As discussed in the previous Section, the vacuum energy of this non-simply connected cavity will have contributions 
coming from the TEM modes, besides the usual TE and TM modes. The existence of TEM modes can be confirmed 
by obtaining explicitly the transverse vector potential A_l defined in Eq(5] that for this particular geometry reads 

Ai = ^- (22) 

The TEM modes contribution to the energy is given by EqfT2l and is independent of r\ and r 2 , as stated above. 

On the other hand, the TE and TM modes contribution do not have such a simple expression and do depend on r\ 
and r 2 . To obtain this contribution, we may start from the already finite expression given in Eq. fT?] and use Cauchy's 
theorem to convert the sum over A^'s and A^'s into a closed path integral in the complex plane of an appropriate 
function. 

If f(z) is a function with "l"-valued simple poles at z — XkD and z = XkN for all XkD,kN then we may compute 
the sum in A = XkD,kN in Eq. [14] as a Cauchy integral, 



E 



XK l{ 2lXa) =2m f z K 1 i^ f[zl (23) 
I Jr I 



if the contour C encloses all the poles of f(z) in z = XkD,kN and the function K\(2lza) is analytic in the interior of 
the curve. The explicit form of f(z) comes out after observing that in the case of the cavity between two perfectly 
conducting concentric cylinders of radius r\ and r 2l the Dirichlct and Neumann cigenfrequencies -which correspond 
to the TM and TE modes, respectively- are all the solutions of (see [l6j ) 

J„(An)iV n (Ar 2 ) - J n (Xr 2 )N Tl (Xr 1 ) = 0, , , 

J4(An)K(Ar 2 ) - j;(Ar 2 )K(An) = 0, ^ 

with n any integer number. (Notice that if A is a solution then —A is also a solution, but since both correspond to the 
same eigenfunction we may keep only the A > solutions to avoid double counting.) From here it is easy to see that 

f(z) = J- z ln [( J n(zri)N n (zr 2 ) - J n {zr 2 )N n {z ri )) (J^(zn)N^zr 2 ) - J' n (zr 2 )N^z ri ))] (25) 

n 

satisfies the above requirements. 

In order to choose the contour C we observe that K\(z) is singular at z = but is analytic for Re(z) > 0, where goes 
to zero as z -1 / 2 e~ Re ^ for large i?e(z). Moreover, since the contour must enclose the real positive axis beginning in 
Xmin (the minimum of the solutions of Eq. [23] for all n) we choose the contour to be a pizza slice with its vertex at 
z = Xmin /2, angle < cj> < n/2, and centered in the real axis: 

z = + pe~^/ 2 p S (0, L) 

C = hm i z = + L e ie 9 e (-<p/2, +<f>/2) (26) 



L — >oc 



z = 



x 2 

A M I N 



+ pe+^Z 2 p e (£,0) 



Once f(z) and C have been correctly chosen, we have an explicit expression for the TE+TM Casimir energy in 
terms of a double sum and a closed path integral, 

E TE +E TM = _ i g /3£l££) /nW) (27) 

i=l,7l=-O0 C 
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where f n (z) is each term in the sum in f(z) (see Eq. [25]). 

To compute explicitly Eq. [57] we need to truncate the I and n sum according to a given precision, and compute 
numerically the integral in the upper and lower segments of C -which are essentially the same-, since the contribution 
in the arc of radius L goes to zero. The criteria used to truncate the sum is best analyzed in Eq. 1141 where the sum 
is exponentially damped by the Bcssel function when its argument grows. In fact, we may divide in Eq. [14] the sum 
in A as different sums for each n (see Eq. I24p . then we define A m ; n (n) as the minimum of the A's for a given n, and 
then for each n we keep I's such that 2lX m i n (n)a < D. Here D is chosen such that all the thrown away terms in the 
sum are damped by at least an e~ D factor. On the other hand, to truncate the sum in n, we set I — 1 and we choose 
n such that 2X min (n)a < D. This criteria should give a precision of order e~ D to the final result in the sum. In our 
calculations we have taken D = 8 which is enough for our purposes. 

It is worth noticing at this point that if we would have performed directly the sum in Eq. [TJ] over all relevant A's 
instead of using the Cauchy integral approach, then we would have had to study the roots of Eq. [M] for each n and 
keep only those which satisfy 21 Xa < D for each I. Although more difficult, this would have also been a possible 
approach. 

The TE+TM Casimir energy has been numerically computed using Eq. [27jfor different cylinders radii (fx and r%) 
and distance between the plates (a) . We have also computed the Casimir force deriving Eq. [27] with respect to a and 
computing numerically the resulting expression: 

oo - 

F TE + F™ = -i V z 2 (K (2lza) + K 2 (2lza))f n (z). (28) 
, , Jc 

l — 1,71— — OO 

As a check for the numerical TE+TM calculation we have corroborated that its behaviour for small and large a 
corresponds to the expected proximity (Eq. [15]) and exponential (Eq. [T6]) behaviours, respectively. In both cases we 
find, as expected, a convergence to unit in the ratio of the numerical energy and its expected asymptotic behaviour. 
In Fig[2]we show this convergence in the proximity limit, which is the more complex from a numerical point of view, 
since it involves the summation of a large number of modes. 

^te+tm [numeric] 



^PFA 




a 



0.5 1.0 1.5 2.0 n 

FIG. 2: Ratio of the numerically computed TE+TM Casimir energy to the expected proximity behaviour (Eq. I15[) in the 
small-a region. Our calculations reach a/ri ~ 0.1 which is enough for our purposes. The lines correspond, from lower to upper, 
to r 2 /n = 2, 1.8, 1.6, 1.55, 1.5, 1.45 and 1.4. 

In order to explore the TE+TM to TEM transition, we have studied the ranges r-z = 1.1 ry to r% = 40 ri and focused 
on the region of a's where the TE+TM energy (or force) equals the TEM energy (or force). As usual, the numerical 
calculation becomes harder when the surfaces become closer, i.e. when r 2 approaches r%. In Figure [3] we have plotted 
the TE+TM Casimir energy as a function of a together with the TEM Casimir energy for three different r 2 /ri 
configurations. As expected, the TE+TM energy dominates in the small-a region, but the TEM energy dominates 
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when the distance between the plates is larger than the critical distance a c . (Note that given r\ and r 2 there are two 
different critical distances, one for the energy and another for the force. With no risk of confusion we use the notation 
a c for both of them.) 




FIG. 3: TEM (dotted) and TE+TM (solid) contributions to the Casimir energy for different ra/n configurations as a function 
of the distance a (in units of ri) between the plates. From left to right the solid lines correspond to n/ri = 1.2, 2 and 4, and 
the corresponding critical distances are a c = 0.44ri, 0.68ri and 1.41 n, respectively. Observe that the TEM contribution is 
independent of Tilr\. 

We have plotted in Figure 2] the Casimir energy and force critical distance a c for several radii ratios r 2 /ri . As it 
can be seen in the plot, the dependence of a c with r 2 becomes rapidly linear. According to Eq ll9i this means that 
flc = T2 g{fi/r2) approaches its limiting value a c rj J"2<7(0) for r 2 > 2r\. In this case we find a c ~ 0.36 r%, for the 
critical distance associated to the energy. 

As mentioned at the end of Section II, the linear relation between a c and r 2 corresponds to the physical situation 
in which the TE+TM Casimir energy of the non-simply connected cavity approaches that of a simply connected one 
of radius r 2 . Therefore, one should be able to obtain the coefficient g(0) by a comparison of the TEM Casimir energy 
with that of a hollow cylindrical cavity of radius r 2 , which is given by Eq l27l with 

f n (z) = ^-\n{J n {zr 2 )J' n {zr 2 )) (29) 

We have checked that this is indeed the case. The Casimir energy for the hollow cavity, obtained again using Cauchy's 
theorem, is plotted in Fig [5] As expected, this energy interpolates between the proximity result at short distances, 
and the exponential behaviour at long distances. The number g(0) is determined by the value of a/r 2 for which the 
energy of the hollow cavity equals — n/12a. In this way we obtain g(0) w 0.36, that coincides with the slope of the 
linear relation between the critical distance for the energy and r 2 presented in Fig. 21 

IV. CONCLUSIONS 

In this paper we have described a geometry in which the existence of TEM modes induce a long range Casimir 
interaction. In particular, we have shown that the electromagnetic Casimir force between two parallel plates inside a 
non-simply connected cylinder is essentially given by the sum of the Casimir forces for 1 + 1 scalar fields with different 
masses. For TE and TM modes, the masses are given by the eigenvalues of the Laplacian on a z = const section of the 
cylinder, with the appropriate boundary conditions, and are non- vanishing. The opposite happens to TEM modes, 
whose zero-point energy corresponds to a massless field, and this is the reason of the different qualitative behaviour of 
their contribution to the force. On the one hand, the TEM force scales as 1/a 2 at all distances and does not depend 
on the area of the plates. On the other hand, in the short distance limit TE and TM forces reproduce the parallel 
plates result proportional to A/a 4 , and in the long distance limit they are exponentially suppressed due to the finite 
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FIG. 4: Critical distance (in n units) where the Casimir energy (lower plot) or force (upper) of the TE+TM modes equals 
that of the TEM modes for different r%jr\ configurations. As it can be seen, the dependence of a c with T2 becomes rapidly 
linear. The slopes of the asymptotes are 0.36 and 0.61. 
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FIG. 5: Casimir energy for two plates separated a distance a in a hollow circular cylinder of radius T2- The expected PFA 
and exponential behaviour for small and large a/r2, respectively, is plotted in dotted lines. The dashed lines show that, as 
expected, this energy equals E TEM = — -n/V2a at ajri ~ 0.36 (see text and Fig. [4|. 



size of the plates (or, in the equivalent picture, to the non- vanishing effective masses). As a consequence, TE and 
TM modes dominate at short distances, while TEM modes do it at long distances. The critical distance where both 
contributions are balanced depends of course on the form of the section of the cylinder, and decreases with its area, 
as we explicitly showed in the particular example of a coaxial circular cylindrical cavity. The summation over the 
effective masses to compute the TE and TM contributions of the force for this particular case was performed using 
Cauchy's theorem, starting from the renormalized Casimir energy for a single massive field in 1 + 1 dimensions. 
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Throughout the paper we considered the Casimir energy for perfect conductors at zero temperature. It would 
be interesting to generalize these results to take into account the combined effects of finite conductivity at nonzero 
temperature. In this context, it is worth remarking that the dominance of the TEM modes over the TE and TM 
contributions to the Casimir force is also valid at a nonzero temperature T, at least for perfect conductors. In fact, 
it has been shown [3] that at sufficiently long distances (Xminci 3> 1 and aT ^ 1), the Casimir force for TE and TM 
modes is proportional to \MiNTe~ 2XMINa , i.e. is exponentially suppressed. On the other hand, in the same situation 
the TEM force is proportional to T/a 

While the existence of the long-range TEM Casimir force is of conceptual interest, it would be very difficult to 
measure it. Indeed, as discussed at the end of Section II, this contribution to the force does not depend on the area 
of the plates, its absolute value is extremely small, and therefore it would be measurable only at very short distances 
with the present technology. However, in this regime the TE and TM forces would be much larger than TEM force, 
unless the area of the pistons is sufficiently small. Therefore, the measurement of the TEM Casimir force should 
involve thin rings at short distances. A rather difficult experiment, indeed. 
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